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Abstract 

In recent work, we used pseudo-differential theory to establish conditions that the initial- 
boundary value problem for second order systems of wave equations be strongly well-posed 
in a generalized sense. The applications included the harmonic version of the Einstein 
equations. Here we show that these results can also be obtained via standard energy 
estimates, thus establishing strong well-posedness of the harmonic Einstein problem in the 
classical sense. 

PACS numbers: 04.20.Ex, 04.25. Dm, 04.25. Nx 
1. Introduction 

In the harmonic description of general relativity, the Einstein equations reduce to a 
constrained system of 10 quasilinear wave equations for the components of the spacetime 
metric. Recently [1] we used the theory of pseudo-differential operators to prove that one 
can construct constraint preserving boundary conditions of Sommerfeld type such that the 
resulting initial-boundary value problem (IBVP) is well-posed in the generalized sense. We 
show in this paper that the decisive estimate can also be obtained by integration by parts. 
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It follows that the full quasilinear system can be treated by standard energy estimates to 
establish that the harmonic IBVP is strongly well-posed in the classical sense [2] . 

Our results have broad application to other systems of second order wave equations 
besides general relativity, e.g. to elasticity theory, acoustics and electromagnetic theory. 
Most analytic and computational treatments of the IBVP utilize the well-developed theory 
of first order symmetric hyperbolic systems. We develop our results here in their natural 
second order form, which avoids the integrability constraints associated with the extra 
variables introduced in a reduction to first order. 

In view of the wide range of potential applications, instead of the geometrical notation 
of general relativity, we present our main results in See's 2 - 3 in a style familiar to a 
broad audience of computational mathematicians and physicists. We use the notation and 
definitions of the classic treatise [2] on the IBVP. In these sections, we treat systems of 
wave equations with constant coeflBcients, as arise in the frozen coefficient form of the 
harmonic Einstein equations. In Appendix 1, we show that our results extend locally in 
time to the well-posedness of quasilinear problems, such as the harmonic IBVP. 

There is an intimate interplay between the treatment of the analytic theory using 
energy estimates and its finite difference approximation [2,3]. The proof of existence of 
analytic solutions outlined in Appendix 1 is based upon a convergent finite difference 
approximation incorporating semi-discrete energy estimates via summation by parts. Con- 
versely, using summation by parts, the energy estimates of the analytic theory can be 
parroted by discrete energy estimates which guarantee the stability of a finite difference 
approximation [3] . In previous treatments of the second order wave equation, the discrete 
energy approach has been used to develop stable difference algorithms for Neumann and 
Dirichlet boundary conditions [4,5] and this treatment has been extended to the Einstein 
equations [6]. The results presented here provide a guide for applying the discrete energy 
approach to second order systems with a wide range of boundary conditions. These include 
the Sommerfeld condition, which has important application to outgoing wave problems, 
but also more complicated conditions involving derivatives tangential to the boundary. 

The set of constraint-preserving boundary conditions allowed for the Einstein system 
is quite extensive. The selection of a "preferred" choice, e.g. Dirichlet, Neumann or Som- 
merfeld, rests upon additional geometrical or physical criteria depending upon the nature 
of the problem. An important example is the description of an isolated radiating system, 
for which the physically appropriate boundary condition should be adapted to the absence 



of external influences. Unless the treatment of such a system is extended to infinity, this 
requires introduction of an artificial boundary for which boundary data may not be known, 
so that it becomes advantageous to use a boundary condition for which homogeneous data 
is a good approximation. We defer this issue to future work where we will present a ge- 
ometric version of the results of this paper to formulate boundary conditions appropriate 
for an isolated gravitational system. Additional subtleties arise in the treatment of bound- 
aries which are moving, such as an oscillating conducting boundary in an electromagnetic 
wave problem or the artificial boundaries that arise in the dynamically curved spacetime 
of general relativity. In Appendix 2, we illustrate for the wave equation on a curved space 
background how the geometric approach can be used to simplify the formulation of energy 
estimates for such moving boundaries. 

2. The main estimate 

Consider the wave equation 

Utt = Ua:x + Uyy + U;^^ + F (1) 

on the half-space 

X > 0, — oo < y < oo, — oo < z < oo, 
with boundary conditions 

aut — Ux + PiUy + P2UZ + oiq, for a; = 0, o: > 0, 01+0^ < o? 1 (2) 

with boundary data q and initial data 

w = /i, ut = f2, t = (3) 

of compact support. The subscripts {t,x,y,z) denote partial derivatives, e.g ut = We 
assume that all coefficients, the data and the solution are real and that a > 0, Pj are 
constants. Also, we use the notation 

{u,v), \\u\\^ = {u,u); {u,v)b, \\u\\% = {u,u)b, 

to denote the L2-scalar product and norm over the half-space and boundary space, respec- 
tively. 
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Strong well-posedness of an IBVP extends the requirements of a well-posed Cauchy 
problem to include estimates of boundary values [2]. In order to adapt the standard 
definition to second order systems, we write u = (u, ut, Ux, Uy,Uz) to represent the solution 
u and its derivatives; and similarly for the initial data we write fi = {fi, fix, fiy, fiz)- For 
the problem (l)-(3), strong well-posedness requires the existence of a solution satisfying 
the estimate 

||u(t)f + \\n{s)\\lds < Kt (llfif + \\f2f + £ \\F{s)fds + Ms)\\lds^ , 

where for every finite time interval <t <T the constant Kt is independent of F, fi, /2 
and q. 

In [1] we used pseudo-differential theory to prove that the problem (l)-(3) is well- 
posed in the generalized sense and that it is boundary stable. The proof in [1] is only 
given in 2D. In 3D the necessary and sufficient condition for a is o; > ^/P^ + This 
same inequality governs the strong well-posedness of the problem. (It corresponds to the 
timelike property of the vector field in the geometric treatment of Appendix 2.) 

We now want to prove that results can also be obtained in terms of standard estimates 
using integration by parts. For simplicity of presentation, we restrict our attention here 
to estimates of the derivatives of u. An estimate for u itself can easily be obtained by the 
change of variable u — > e^^u, as described in Appendix 1. We start with 

Lemma 1. Let 71, 72 be real and 7^ -|- 7! < 1. Then 

E = \\ut\\^ + \\ux\\^ + ||%||^ + \\uz\\^ - 2(ttt,7i% -I- 72W«) (4) 
is a norm for the derivatives (ut, Ux,Uy,Uz). 

Proof. Since 7^ -|- 7! < 1, we can choose 5 with < 5 < 1 such that 

2|(wt,7i% + 72^2)1 < (1 - 5)||wtf + Y^hi +72)i\Wyf + hzf) 

<{l-6) {\\Utf+\\Uyf+\Kf). 

This proves the lemma. 

Let 7i = 72 = Now we can prove that there is a standard energy estimate. 
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Theorem 1. The solution of (l)-(3) satisGes the energy estimate 

+ <E+ \\Ff + a\\q\\l. (5) 

Proof. Integration by parts gives us 
d d 

Q^htf = 2{ut,utt) = -Qlihxf + \\uyf + IM'^) + 2{ut,F) - 2{ut,Ua,)B (6) 
d 

where, for example, we have used {uyy,Uz) = —{uy,Uyz) = 0. Since (2) imphes 

2 

2{ut,Ua;)B = -\\UxfB + 2{Ua;,-fiUy + -f2Uz)B + 2{Ua;,q)B, 

by subtracting (7) from (6) we obtain 

d 2 

—E = 2{ut - 71% - -i2Uz,F) - -||"a;||| " 2{Ux,q)B 

< \\ut - 7i% - 72%|P + ll-P'lP - + 



The identity 



\\ut - 7l% - l2Uz\\^ = E - \\Ux\\^ - \\Uy\\^ - \\Uz\\^ + ||7l% + l2Uz\\^ 

then imphes (5) and thus proves the theorem. 
The theorem teUs us that we can estimate 

£;(T) and f \\ux\Wdt in terms of E{{)), [ \\F\fdt and f \\q\\%dt. 
Jo Jo Jo 

For the apphcation to the Einstein equations we also need estimates of the boundary norms 
of Uy and Uz- We have 

d 

— (Wa;, Ut) = {U^t, Ut) + (Wa;, 

= -^ll^lll + {Ux.U:,^) + {u^^Uyy) + {U:,,Uzz) + {u^^^F) (8) 
= -^Iktlll - \\\Ux\\% + \\\Uy\\B + \\\Uz\\b + {U^^F). 
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The boundary conditions (2) give us, for any 6 with < S < 1, 

ll^tlll = Il7i% + l2Uz + -Ux + q\\% 

a 

1 1 

< WllUy + 72W2||| + '^WllUy + 'y2Uz\\B\\-Ux + q\\B + \\—^x + qWe 

< (1 + S)\\^iUy + 72«^||| + (1 + ^) II ^Ux + q\\l 

<{l + 6) (7f + 7l) + ll^.lll) + (1 + bW-Ux + q\\l. 

da 

Since ^ := 7i + 7I < 1 we can choose 5 such that {1 + 6)g < {1 — 5). Therefore, by (8), 

11 d 
^ihyWl + hzWl) < (1 + ■t)\\-Ux + q\\l + hxWl + 2—{ux,Ut) - 2{ux,F). 

da at 

Since {uxi Ut) can be estimated by E, we have proved 
Theorem 2. 



/ (ll'^tlll + \M\% + \\Uy\\B + \\'^z\\%) dt 
Jo 



T pT 

2j+ , / ll„l|2 



< const. (£;(0) + / \\F\\^dt+ / \\q\\%dt). 
Jo Jo 

The results can easily be generalized to half-plane problems for wave equations of the 
general form 

utt — PoUt + Piu, xi > 0, —oc < Xj < 00, j — 2,3 (9) 
with boundary conditions 

aUt = Wa;i + PlUx^ + ^2^X3 , Xi =0. (10) 

Here 

3 3 
Pq — Cjd/dxj, Pi — Gijd^ /dxidxj, an > 

j = l i,j = l 

and Pi is strongly elliptic. 

By general coordinate transformations, we can transform (9), (10) to the simple prob- 
lem (l)-(3), except if the coordinate system moves normal to the boundary. To discuss this 
case we consider a boundary moving with constant velocity c/y/l + in the x-direction 
and transform (1) according to 

X |~ d III 

X = r- ^ , t = y = y, z^z, 

vi + 
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so that the moving boundary is located at x' = 0. After dropping the primes, we obtain 



and we consider the problem (l'),(2),(3). The same pseudo-differential technique as in 
Section 2 of [1] shows that this new problem is well-posed in the generalized sense if and 
only if 



Subject to this inequality, in Appendix 2 we use a geometric approach to establish strong 
well-posedness for this moving boundary problem by the energy method. 

3. The Einstein equations 

We consider the half-plane problem treated in [1] for the harmonic Einstein equations, 
which we know describe. In the frozen coefficient formalism based upon an orthonormal 
frame, the components of the densitized spacetime metric satisfy the wave equations 



Uu = 2cUxt + Uxx + Uyy + Uzz+ F 



(1') 




( 




= F 



(11) 



X > 0, t > 0, — oo < y < oo, — oo < ^ < oo. 



where F consists of lower order terms. 



These equations are subject to the constraints 




(l2) 



Because the constraints satisfy homogeneous wave equations, if C°'{0,x,y, z) = and 
^C"*(0, X, y, z) = 0, a = (t, x, y, z), then they remain zero at later times if = are part 
of the boundary conditions for (11) at a; = 0. 
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As in [1] , we consider the choice of boundary conditions 
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where q represents the boundary data. Here ai, 02, 03, hi, 62 and ci, C2 are real constants 
such that the eigenvalues Xj of 



(14) 



are real and negative. The constraints vanish on the boundary provided g*" = 0, a = 
{t,x,y,z). 

We can now obtain standard energy estimates for the Einstein equations both in the 
interior and on the boundary. The following theorem strengthens the results in Section 3 
of [1] to strong well-posedness: 
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Theorem 3. The half-plane problem for the system (11) with constraints (12) and bound- 
ary conditions (13) is strongly well-posed if the eigenvalues of the matrices (14) are real 
and negative. 

Proof. As explained in Appendix 1, it suffices to show that all the first derivatives are 
bounded. We start with the last components. Since the boundary conditions do not contain 
terms in the tangential direction, Theorems 1 and 2 tell us that all the first derivatives of 
7^^ can be bounded in terms of ||F|| and \\q\\B- Thus we gain one derivative both on the 
boundary and in the interior. The same result holds for the derivatives of 7^^ and 7^^. 

In the same way as in [1], the boundary conditions for 7*^, 7^^ can be decoupled by a 
unitary matrix U such that 

Introducing new variables by 




we obtain the equations 

dt\ry \o X2jdx\r'j W"J 

where 

(?)=K-I^Co>ICo> (:::))■ 

We have already estimates of q^^ and g^^^. Therefore we can estimate all the first derivatives 
of 7^^ in terms of ||F|| and \\q\\B- The same is true for 7*^. This process can be continued 
for the remaining components. Thus we have proved Theorem 3 of [1], where we can 
remove "in the generalized sense" . 

Theorem 3 is also valid when the matrices of (13) for the tangential derivatives are 
upper triangular, i.e. only terms above the diagonal are not zero (or equivalent to that 
form by unitary transformation). This allows the sequential argument in the proof. It also 
generalizes to full matrices which are sufficiently close to upper triangular form. A fuller 
discussion of the most general case will be given in future work. 

Local existence theorems for quasilinear equations follow by iteration of the linearized 
equations, as described in Appendix 1. These results establish strong well-posedness, 
locally in time, of the quasilinear harmonic IBVP. 
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4. Discussion 

We have shown that the results of [1], obtained using pseudo-differential theory, can 
also be derived in a more transparent way based upon integration by parts to establish 
strong well-posedness of a broad class of IB VP's governed by second order quasilinear 
wave equations. The underlying arguments require no need to rewrite the wave equations 
as a first order system. As shown in Appendix 1, for smooth data there exist estimates 
for arbitrarily high derivatives (which is a key requirement for treating the quasilinear 
case). The boundary conditions are flexible and are stable against perturbation of their 
coefficients. These properties are important for numerical calculations. 

We anticipate that these results will have application to the broad class of problems 
based upon second order wave equations. For the standard wave equation (1), the chief 
restriction is that the coefficient a of ut in the boundary condition (2) is positive and 
satisfies > 0^+0^; for the general wave equation on a curved space background, 
the corresponding restriction is that the vector T" introduced in Appendix 2 is timelike 
and future directed. This encompasses boundary conditions of Sommerfeld type. The 
application to Maxwell's equations expressed in terms of a vector potential in the Lorentz 
gauge is straightforward. 

A prime motivation for this work is the formulation of a well-posed IBVP for the 
system of harmonic Einstein equations. The importance of numerical simulations in general 
relativity has spurred a large number of works which have established many of the necessary 
ingredients for a well-posed IBVP. For a review see [7]. The first complete well-posed 
formulation was given by Friedrich and Nagy [8] for a version of Einstein's equations 
in which the curvature tensor, i.e. quantities constructed out of second derivatives of 
the metric, was included in the evolved variables. Choquet-Bruhat's first proof of the 
well-posedness of the Cauchy problem for Einstein's equation was given in the harmonic 
formulation [9]. It is satisfying to extend her work to the harmonic IBVP. Our results 
also immediately apply to the generalized harmonic formulation in which harmonic forcing 
terms are allowed [10], since this does not change the principle part of the system. 

Since the pioneering results of Pretorius [11,12], there has been rapid progress in the 
development of numerical codes based upon the generalized harmonic formulation with 
the capabilty of simulating relativistic binaries consisting of black holes or relativistic stars 
[13-15]. This difficult problem requires extensive computational tools that do not enter 
the analytic treatment considered here, namely grid refinement and numerical dissipation. 
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Nevertheless, our results provide the analytic background necessary to justify a successful 
numerical treatment and, at least for finite difference codes, they provide further guidance 
on how to establish robust boundary conditions. Examples of the boundary conditions 
considered here have been incorporated in a unigrid, second differential order, harmonic 
code and successfully tested on model problems [16]. It is beyond the scope of the present 
work to suggest how they might be incorporated in specific black hole codes. However, 
in future work we will show how the geometric approach of Appendix 2 can be further 
developed to yield boundary conditions which are well tailored to the treatment of isolated 
astrophysical systems. 

Appendix 1. The quasilinear case 

The preceding energy estimates establish that a solution of the IBVP with frozen 
coefficients is unique and depends continuously on the data. In this section we want to 
show that local existence theorems and energy estimates for second order quasilinear wave 
equations are proved in the same way as for first order symmetric hyperbolic systems. It 
all depends on a priori estimates for arbitrarily high derivatives of the solutions of linear 
equations with variable coefficients. Consider the halfplane problem for 

Utt = Pu + Ru + F, X > 0, —oo < y < oo, (15) 

with boundary conditions 

aut = Ux — + ru + q, a{y, z,t) > S > 0, (16) 

and initial data 

u{t^O)^h, ut{t = 0)^/2. (17) 

Here 

Pu = {aux)x + {buy)y - 2^ut - 7^tt 

and 

Ru = CiUt + C2Ux + C^Uy + c^u. 

Ru are terms of lower (first and zeroth) differential order. All coefficients are smooth 
functions of a;, t and a > ao > 0, 6 > 6o > 0, ao, 6o and 5 are strictly positive constants. 
The initial data are smooth functions which are compatible with the boundary conditions. 
Here 7 > is a constant obtained by the change of variables u — > e^^u' and then deleting 
the 'prime'. This introduces the term 7^||tt|p in the energy E in (18), which provides an 
estimate of llttlP. 
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Lemma 2. There is an energy estimate which is stable against lower order perturbations. 



Proof. Integration by parts gives 
dt ' dt 



d d 



(18) 



= -47||'Ut||2 + 2{ut, F) + 2{ut, Ru) - 2{ut, au^)B 

+ at||Wa;||^ + 

< const. + - 2(ut,aUa;)B. 

Using the boundary conditions gives 

-{ut, aUx)B = -{ut, aaut)B - {ut, 'yau)B + (ut, ru + q)B 

= -{ut, aautjB - {ut, ^aou)B - {ut, 7(0 - ao)u) ^ + {ut, aru + aq)B 

Id 1 
< -2^cio-Q^\H\B - 2 ^(^'^t)B + const. {\\u\\% + WqWb)- 

Therefore, we obtain from (18) 

d 

— {E + 'yao\\u\\% + {ut,aaut)B) ^^g^ 
< const. (-E + ||tt||| + + ||g|||) . 

This proves the lemma. 

Now we can estimate the derivatives. Let v — Uy,w — Uf. Differentiating the differen- 
tial equation gives us 

Vtt = PV + RV + RyU + {ayUx)x + {hyV)y + Fy, 

(20) 

wtt = Pw + Rw + Rtu + {atUx)x + {hv)y + -Ft, 

respectively. 

RyU and Rtu are linear combinations of first derivatives of u which we have already 
estimated and can be considered part of the forcing. 
The differential equation (15) tells us that 

o-Uxx = wt — bvy + terms we have already estimated. 

Thus Uxx is lower order with respect to v, w and, except for lower order terms, v, w are 
solutions of the same differential equation as u. The same is true for the boundary con- 
ditions. Therefore we can estimate all second derivatives. Repeating the process, we can 
estimate any number of derivatives. 
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We can now proceed in the same way as in [2], where we have considered first or- 
der systems, to obtain existence theorems for equations with variable coefficients. We 
approximate the differential equation by a stable difference approximation and prove, us- 
ing summation by parts, that the corresponding estimates for the divided differences hold 
independently of the gridsize. In the limit of vanishing gridsize, we obtain the existence 
theorem. Since we can estimate any number of derivatives, it is well known, using Sobolev's 
theorem, that we can obtain similar, although local in time, estimates for quasilinear sys- 
tems. By the same iterative methods as for first order symmetric hyperbolic systems, it 
follows that well-posedness extends locally in time to the quasilinear well as other 

standard results such as the principle of finite speed of propagation. 

Appendix 2. Geometric Derivation of the Estimates 

We now use a geometric approach to show how the estimates established in Sec. 2 
can be extended to the general initial-boundary value problem for the second order wave 
equation on a curved spacetime. A similar geometric approach has been used to implement 
the second order harmonic formulation of Einstein's equations as an evolution-boundary 
code based upon summation by parts [16]. See also [17] for another application of the 
geometric approach to the treatment of boundaries. 

Using standard notation of general relativity, we consider the wave equation 



for a massless scalar field propagating on a Lorentzian manifold with boundary of the form 
M = [0,T] X S, where S is a compact, 3-dimensional manifold with smooth boundary 
dT,, each time-slice = {t} x E is spacelike and the boundary T = [0,T] x dT, is 
timelike. Here V denotes the covariant derivative associated with the spacetime metric 
g with signature (— , -|-, -h, -h). The initial-boundary value problem consists in finding 
solutions of (21) subject to the initial conditions 



^""VaVb^ = F 



(21) 



with Cauchy data / and h ohHq, and the boundary condition 



[{T^ + aN^)Wb(l>] 



T 



= Q 
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with data q on T. Here n and N° denote the future-directed unit vector field to the time- 
shces T^t and the outward unit normal vector field to T, respectively; is an arbitrary 
future-directed timelike vector field which is tangent to the boundary surface T; and a > 0. 
The motion of the boundary is described geometrically by the hyperbolic angle N^rn,. 



— 1 . A Sommerfeld boundary condition then corresponds to the choice a = 1 for which T'' -|- 
points in an outgoing characteristic (null) direction picked out by the geometry. See [16] 
for a numerical study of the harmonic Einstein system carried out with such a Sommerfeld 
condition, with evolution in the direction. In [16], the choice of matrices governing 
the tangential derivatives in the boundary conditions (13) was made for mathematical 
simplicity. In future work, we will explore more physically and geometrically motivated 
choices. 

In order to establish estimates, we introduce the notation (f)a = ^a(f> and the energy 
momentum tensor of the scalar field 



The essential idea is the use of an energy associated with a timelike vector u"' = T"' + 5N°', 
where (5 > 0, so that it" points outward from the boundary. The corresponding energy 
E{t) and the energy flux J^{t) through the boundary are 



Without loss of generality, we can assume that is normalized such that gbcT^T^ 



©6 = 06<^' 



a 




and 




It follows from the timelike property of that E{t) is a norm for (pait)- 
Energy conservation for the scalar field, i.e. integration by parts, gives 




E 



SO that 




(22) 



E 
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The required estimates arise from considering the flux density 

where ffbc = 9hc + Tfo^c — Ni,Nc is the positive definite metric in the tangent space of the 
boundary orthogonal to T**. By using the boundary condition to ehminate T°'(f)a in the 
last group of terms, we obtain 



SO that 



+ {-a + 6{1 + a^)) (iV>j2 ^ _ 2aS)N''(l>aq + 5q^ 



+ (-a(l-e) + 5(l + a2))(iV-(/,,y 
-ea N''(j)a - ^ -q + U + ^ ' 



2ae J \ 4ae 
In the above equation, we have introduced the e-terms, with < e < 1, in order to establish 
the inequality 

+ (-a(l - 6) + 5{1 + a^)) {N-cl>af (23) 



4ae 

The boundary estimate of 00 now follows by requiring — a(l — e) + 5(H-a^) < 0, which 
guarantees that is timelike. With the choice 

a(l — e) 



(22) and (23) give 



dtE^ J ^(^(Ar>,)2 + (T>o)2 + iy«Va0i 



(24) 



< const. {E+ I + I q^ 



E as 
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In analogy with (19), this leads to the required estimate of the gradient 0a on the boundary 
(as well as the usual estimate of the gradient at a fixed time) to prove that the problem 
is strongly well-posed. An estimate of itself follows by introducing a mass term in (21) 
through the change of variable — > e'''*^', as described in Appendix 1. Energy estimates 
for the problem (l'),(2),(3) in Sec. 2 follow from the choice T*' = {T^T^^ ,T^) = 
t{q' , 0, — /S^, —132) with = —'ug^°'VaX = —i^{c, 1, 0, 0), where r and v are positive nor- 
malization constants and a' and (31 are related to the coefficients a and (3i in (2) by 

ra' 

a = c 

av 

Pi^^, i = l,2. 
au 
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